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ABSTRACT 

Simultaneous  measurements  of  waves,  orthogonal  water 
particle  velocities,  temperature,  conductivity  and  acoustic 
phase  and  amplitude  measurements  at  20  and  40  kHz  over 
horizontal  and  vertical,  2m  paths  were  made  at  the  NUC 
oceanographic  tower,  San  Diego  in  April  1974. 

Four  stochastic  models  for  acoustic  propagation  including 
geometric,  Born,  Rytov  and  Debye  approximation  models  were 
tested.   The  Krosil'nikov  -  Obhukov  formulation  of  the  Born 
approximation  yielded  reasonable  results.   Difficulty  in 
establishing  scales  obviated  good  results  from  the  Rytov 
and  Debye  models. 

Taylor  thermal  microscales  from  spatial  correlation 
analyses  determined  vertically  across  the  thermocline  varied 

from  6.8  to  178  cm.   The  measured  variance  of  the  index  of 

-4  -3 

refraction  ranged  from  4.5x10   to  2.7x10   . 

Temporal  stochastic  analyses  showed  horizontal  acoustic 

path  losses  were  ten  times  vertical  losses  because  of 

diffraction  effects.   The  primary  mechanisms  for  acoustic 

variability  were  internal  and  surface  waves  pumping  the 

thermocline  up  and  down. 
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I.       INTRODUCTION 

A.  HISTORY 

A  continuing  program  of  in  situ  ocean  microstructure 
research  has  been  conducted  at  the  United  States  Naval 
Postgraduate  School  (NPS)  since  1971.   Under  this  program, 
a  series  of  three  annual  experiments  dedicated  to  the  study 
of  the  complex  temporal  and  spatial  relationships  between 
multiple  parameters  measured  in  shallow  water  have  been 
conducted  by  the  Department  of  Oceanography.   Each  of  these 
experiments  has  been  restricted  in  scope  to  measurements 
inside  a  9  mJ  volume  of  water  and  to  oceanographic  phenomena 
which  occur  in  the  frequency  range  0.001  to  2.5  Hz. 

B.  OBJECTIVES 

In  this  thesis,  initial  analysis  of  certain  parameters 
measured  in  the  latest  of  these  experiments  is  reported. 
Primary  emphasis  is  given  to  stochastic  temperature  relation- 
ships affecting  high  frequency  (20  kHz)  CW  sound  transmission 
over  a  short  path.   The  particular  objectives  of  this 
experiment  were  to  investigate  the  dependence  of  acoustic 
phase  and  amplitude  fluctuations  on  the  magnitude  and 
characteristics  of  the  ocean  thermal  structure  and  the 
intensity  of  the  related  turbulence. 
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II.   THEORY 

The  description  of  oceanic  turbulence  by  stochastic 
techniques  has  followed  developments  in  the  atmosphere. 
One  of  the  principal  difficulties  in  application  of  the 
theoretical  methods  provided  from  atmospheric  research  has 
been  determination  of  the  range  of  applicability  of  certain 
turbulence  models  in  terms  of  ocean  parameters.   Another 
has  been  the  application  of  these  models  to  sound  propagation 
in  a  manner  similar  to  the  relationships  developed  for 
electromagnetic  wave  propagation  in  the  atmosphere  as 
summarized  by  Strohbehn  [1968].   Historically,  attempts  to 
resolve  these  difficulties  have  been  interconnected  since 
their  primary  relevance  is  to  acoustic  prediction  which 
motivated  the  research. 

A.   STOCHASTIC  MODELS  FOR  THE  MEDIUM 

Three  expected  value  techniques  which  have  been  applied 
to  turbulence  calculations  in  the  ocean  are  1)  the  spatial 
correlation  (autocorrelation)  function,   2)  the  structure 
function  and   3)  a  "piecewise  correlation"  function. 
Application  of  each  technique  requires  a  statistically 
homogeneous  medium  but  the  order  above  indicates  their 
relative  applicability  to  increasingly  less  stringent 
stationarity  requirements.   The  function  form  for  each  is 
listed  below. 
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The   autocovariance    function,    R   (p): 


Rx(p)    =   E[T(x)T(x+p)]      .  (1) 


The  structure  function,  D  (p) 


Dx(p)  =  E{[T(x)  -  T(x+p)]2}   .  (2) 


The  piecewise  correlation  function,  R'(p) 

A. 

R'(o5  =  ET1  -  E([T(x)  -  T(x+p')]2> 

*Vp;    hL1 2  R  (0) 

x 


D  (pM 

- E^  -  2-fror]     •  <3) 

x 


E  indicates  expected  value  and  p  is  the  spatial  lag. 

For  completely  stationary  processes,  the  structure  function 

and  the  "piecewise  correlation"  function  reduce  to 


Dx(p)  =  2[Rx(0)  -  Rx(p)l   ,  (*») 


and         Rx(p)  =  E[l  -  Dx(  p)  V2RX(0)  ]  =  Rx(pf)/Rx(0) 


(5) 


The  stationarity  difficulties  often  encountered  in  the 
use  of  R  (p)  are  partly  overcome  by  the  use  of  D  (p)  because 
the  large  scale  effects  in  the  ocean  are  differentially 
filtered  out  [Black  1965].   The  piecewise  correlation  function 
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can  be  applied  to  measurements  which  have  only  "piecewise 
stationarity. "   This  makes  it  the  most  powerful  tool  of  the 
three  techniques.   "Piecewise  stationarity"  is  defined  as 
stationarity  (in  the  ordinary  sense)  over  a  sufficient 
record  length  for  which  the  structure  function  and  variance 
can  be  calculated.   This  method  was  proposed  by  Whitmarsh, 
Skudrzyk  and  Urick  [1957]. 

Various  analytical  equations  have  been  used  to  approxi- 
mate the  form  of  the  auto  covariance  function  used  to  describe 
the  environment.   Liebermann  [1951]  used  a  simple  exponential 
form  in  his  early  experiments.   Others  have  used  a  Gaussian 
form  with  some  success  [Stone  and  Mintzer  1953].   Use  of  the 
exponential  form  yields  the  best  fit  to  experimental  data 
above  certain  spatial  lags,  but  it  underpredicts  the  number 
of  "small"  turbulent  blobs.   This  undesirable  effect  is 
eliminated  by  the  use  of  the  Gaussian  form. 

Kolmogorov  [19*11]  predicted  a  2/3  power  law  dependence 
for  the  structure  function  for  turbulent  flow  in  which  there 
was  a  decoupling  of  the  energy  producing  and  energy  dissipa- 
tion regions.   Batchelor  [1953]  transformed  this  result  and 
demonstrated  that  in  wave  number  space  energy  spectra  have 
a  wave  number  to  the  -5/3  power  law  dependence.   This  was 
first  verified  convincingly  by  Grant  et  al .  [1962]  for  ocean 
measurements  and  Pond  et  al.  [1963]  for  atmospheric  measure- 
ments.  Ocean  measurements  do  not  always  show  the  dependence 
since  the  ocean  tends  to  be  stratified  much  more  intricately 
structured  than  the  atmosphere  which  makes  the  recording  of 
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continuous  stationary  data  unlikely.   However  the  -5/3 
power  law  depends  not  so  much  upon  isotropy  as  upon  the 
clear  separation  of  scales  of  energy  production  and  energy 
dissipation  in  which  a  cascade  of  energy  from  larger  to 
smaller  scales  occurs.   Whether  or  not  such  a  form  is 
ubiquitous  in  the  ocean  is  not  as  significant  as  the  presence 
of  horizontal  sheets  which  could  cause  multiple  reflection 
layers.   Very  sharp,  multiple  thermal  layers  have  been 
reported  in  the  Arctic  [Denner  1971],  for  example. 

Recently  it  has  been  shown  that  surface  wave  effects  can 
be  removed  from  ocean  turbulence  measurements  with  the  result 
that  purely  turbulent  flow  remains  and  a  -5/3  power  relation- 
ship may  be  observed  [Thornton  and  Boston  197*0.   These 
findings  were  based  upon  measurements  both  in  and  around  the 
thermocline  and  showed  that  anisotropy  was  not  a  critical 
difficulty,  provided  additional  energy  sinks  or  sources  did 
not  result  (e.g.  in  the  thermocline). 

Unfortunately,  other  energy  contributions  to  the 
turbulence  spectra  do  exist.   Microstructure  related 
turbulence  from  salt-fingering  is  one  form  [Gregg  1973 , 
Denner  1971].   Billow  turbulence  resulting  from  shear 
instabilities  related  to  internal  waves  has  also  been 
reported  by  a  number  of  investigators  [Woods  and  Wiley  1972, 
Gregg  1973,  Powell  197^].   Further,  the  existence  of  the 
internal  wave  itself  could  be  expected  to  add  energy  at  low 
frequencies  as  could  surface  wave  interaction  with  the 
thermocline . 
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B.   WAVE  PROPAGATION  MODELS 

Four  methods  derived  from  the  Helmholtz  equation  have 
been  used  to  describe  wave  propagation  through  the  turbulent 
ocean.   These  are:   1)  geometrical  optics,   2)  the  Born 
approximation  [Born  1933],   3)  Rytov's  method  [Rytov  1937], 
and  4)  the  Bebye  approximation  [Frisch  1968].   While  it  is 
not  a  necessary  assumption,  most  applications  of  these 
methods  presume  that  the  ocean  is  homogeneous  so  that  the 
scalar  wave  equation  can  be  used.   Such  an  assumption  places 
a  rather  severe  restriction  on  the  acoustic  wavelength,  A, 
such  that 


X  <<  zq  (6) 


where  it  =  inner  scale  of  the  turbulence.   This  restriction 
o 

will  be  discussed  in  detail  later. 
1.   Geometrical  Optics 

The  geometrical  optics  (ray  theory)  approach  is  a 
linearization  of  the  wave  equation  which  assumes  only 
Rayleigh  scattering  and  refraction  effects  are  important. 
Since  diffraction  effects  are  ignored,  the  high  frequency 
(short  wavelength)  restriction, 


X  <<  lQ2/L  (7) 


where  L  =  acoustic  path  length,  is  placed  on  the  results. 
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One  of  the  first  experiments  which  used  a  stochastic 
technique  to  examine  the  acoustical  effects  of  small-scale 
temperature  inhomogeneities  in  the  ocean  was  conducted  by 
Liebermann  [1951].   Through  the  use  of  the  classical  optical 
model,  Liebermann  was  able  to  predict  certain  scintillations 
in  sound  amplitude  observed  at  sea  by  the  equation: 

vi2  =  T5  **  °-2     (L/a)3  (8) 


VI2  =  E^>2]  (9) 


a2   =  E[(^-)2]  (10) 

co 

2 
where  Vy   =  variance  of  the  sound  intensity,  I  =  mean  sound 

intensity,  AI  =  the  instantaneous  deviation  of  sound  intensity 

p 

from  the  mean,  a   =  variance  of  the  index  of  refraction, 

L  =  acoustic  path  length  and  a  =  Taylor  temperature  micro- 
scale  [Lumley  and  Panofsky  1964].   Liebermann  also  showed 
that  sound  reverberation  encountered  in  the  ocean  was  larger 
than  that  predicted  due  to  thermal  inhomogeneities  alone. 
This  was  substantiated  by  later  developments  which  have 
shown  that  biological  scattering,  chemical  relaxation  and 
bubbles  in  the  near-surface  region  are  important  contributors 
to  sound  attenuation. 

Sagar  [1955,1957,1959]  was  able  to  show  that  sonar 
amplitude  fluctuations  resulting  from  temperature  measure- 
ments could  be  correlated  to  sea  state.   Additionally,  he 
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showed  that  the  linear  dependence  of  the  variance  of  sound 
pressure  level  on  range  [Sheehy  1950]  was  apparent  only  when 
thermal  microstructure  was  present. 

A  modified  geometrical  optics  technique  has  recently 
been  proposed  by  Tatarski  [1967].   The  method  is  similar  to 
the  classical  one  except  the  path  length  is  described  as 
the  eikonal  or  total  path  length,  6.   However,  one  would 
assume  that  this  would  modify  the  restriction  on  the  results 
to 


X   «   I   2/6  (11) 


No  application  of  this  method  to  the  ocean  has  yet  been  made. 
2.   Born  Approximation 

The  Born  approximation  is  a  perturbation  technique 
applied  to  the  wave  equation  which  is  valid  only  as  long 
as  the  acoustic  amplitude  fluctuations  remain  small.   This 
is  not  a  serious  restriction  in  the  ocean,  even  under  the 
most  severe  scattering  conditions,  but  the  model  only  accounts 
for  a  single  scattering  from  a  turbulent  blob  over  the 
acoustic  path  length,  L.   This  requires  that 

X    «    a  /aL  (12) 

The  geometrical  optics  results  are  a  special,  high  frequency 
case  included  in  this  model. 
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Stone  and  Mintzer  [1962,  1965]  conducted  a  series 
of  laboratory  experiments  from  which  they  demonstrated 

solutions  to  the  Born  Model  which  predict  the  variance  of 

p 

normalized  sound  pressure  level,  V  ,  in  a  low  frequency 

region  (wave  region)  and  in  the  ray  theory  region.   For 
the  wave  region  they  showed: 


V2  =  J  /tt  k  2  a2   a  L   (Gaussian  R  (p)   (13) 


x 


V2  =  a2  k   a  L    (Exponential  R  (p)    (1*1) 


V2  =  E[(^)2]  (15) 

o 


X  <<  2ua  (16) 


X  <<  2ttL  (17) 


kQ  =  2tt/X  (18) 


Where  P   =  mean  sound  pressure  level  and  AP  =  pressure  level 
o 

difference  from  the  mean.   For  the  ray  region,  no  solution 
is  possible  if  the  exponential  form  of  the  spatial  correlation 
function  is  used  since  the  solution  is  divergent.   But  they 
showed  that  the  following  equation  results  from  use  of  the 
Gaussian  correlation  function: 
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y2  =  T5    ^   ct2(L/a)3  (19) 

Neubert  and  Lumley  [1970]  used  data  from  the  above 
experiments  to  compare  with  atmospheric  data  based  upon  the 
same  model.   Their  results  showed  that  the  turbulent  scales 
for  air  and  water  are  of  the  same  order  of  magnitude  for 
"similar"  turbulent  flows.   They  were  also  able  to  better 
define  the  limit  on  wave  length  specified  in  relation  (6), 
such  that 

1 
X  <  a/5[ak(u/co)]2  (20) 


a2  =  Xg2/ak  (21) 


Rt2  *  ^(rig/Xg)4  (22) 


where 


X   =  Taylor  turbulence  microscale, 
g 


a.  =  Prandtl  number  for  the  medium. 


u  =  local  turbulent  velocity, 


c   =  mean  speed  of  sound  in  the  medium, 
o 


R  =  Richardson  number, 
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and 


n  =  Kolmogorov  turbulence  microscale, 

D 


In  addition,  they  found  a  restriction  on  acoustic  frequency 
was  necessary,  though  it  was  usually  less  restrictive  than 
(20).   This  restriction  is 

1  1 
f  I   750[ak(Rt)2]2/Xg   (sec-1)  (23) 

It  follows  from  (6)  and  (20)  that 


lQ   z   0{a/50[ak(u/co)]2}    (cm)  (24) 


which  lends  additional  meaning  to  relations  (6),  (7)  and 
(11). 

Whitmarsh,  Skudrzyk  and  Urick  [1957]  conducted 
experiments  in  the  ocean  which  largely  agree  with  Stone's 
and  Mintzer's  laboratory  results.   However,  they  noted  the 
variance  in  sound  pressure  level  due  to  thermal  scattering 
varied  directly  as  the  acoustic  frequency  and  as  the  range 
squared  in  the  Kolmogorov  turbulence  range.   This  result  is 
obviously  in  the  intermediate  range  between  the  ray  and 
wave  solution  regions  suggested  by  Stone  and  Mintzer. 

They  also  observed  that  the  turbulent  patches  were 
not  spherical,  but  were  spheroidal  in  shape.   In  this  case, 
a  simple  transform  to  cylindrical  coordinates  shows  that  the 
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Taylor  microscale,  a,  must  be  replaced  such  that 


a2  =  (b2/3)(2  +  c2/b2)  (25) 


where  b  =  narrower  dimension  and  c  =  longer  dimension  of 
the  turbulent  patch.   This  has  the  effect  of  reducing  the 
focusing  length  of  an  individual  turbulent  "lens"  which, 
in  turn,  increases  the  effective  ray  solution  region. 

The  following  expression  for  the  variance  of  the 

2 
sound  pressure  level,  V  ,  has  been  presented  as  valid  over 

the  entire  wave  and  ray  regions  [Krasil'nikov  and  Obhukov, 

1956[: 


V2    =  kQ2   a2   a  L[l   -   1/Dtan   1D]  (26) 


D  =   TrL/koa.  (27) 


For  this  result,  the  wave  region  is  defined  when 


X   »   2a2/L  (28) 


and  the  ray  region  is  defined  for 


A  <<  2a2/L.  (29) 
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3.   Rytov's  Method 

Rytov's  method  (method  of  smooth  perturbations)  is 
the  method  used  by  Tatarski  [1Q61]  to  develop  a  stochastic 
wave  model  for  the  atmosphere.   It  Includes  the  geometrical 
optics  and  Born  approximation  methods  as  special  cases  and 
includes  multiple  scattering  effects  as  well.   Unfortunately, 
the  scattering  effects  are  included  such  that  they  cannot 
be  parameterized  directly  in  terms  of  the  turbulence. 

Either  the  spatial  correlation  or  structure  function 
may  be  used  with  this  method  as  dictated  by  considerations 
regarding  the  stationarity  of  measurements  to  be  analyzed. 
The  restriction 


«v 


£  Vl  (30) 


is  necessary  for  reasonable  solutions  of  the  resultant 
equations.   It  is  also  necessary  to  compute  the  range  of 
validity  of  the  perturbation  expansion  when  this  method  is 
used.   This  calculation  will  not  be  discussed,  but  Strohbehn 
[1968]  suggested  that  the  method  is  valid  as  long  as  the 
variance  in  the  normalized  log-amplitude  of  the  sound  pressure 
level  is  less  than  0.8.   He  also  argued  that  the  results 
from  sound  phase  calculations  should  be  less  restricted  than 
those  for  amplitude. 

Black  [1965]  applied  this  model  to  certain  data  from 
the  CHAIN  experiment.   One  important  conclusion  he  reached 
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was  that  the  standard  deviation  of  the  fluctuation  in  sound 
intensity  through  a  turbulent  medium,  oa,  and  the  principal 
period,  T  ,  of  its  fluctuation  can  be  predicted  by  the  equations 

11  _  _L 

as   =   12.7   Cn  L12    X     12  (3D 

Ts   =    7r/XE/2vn  (32) 

C2n   =   CT2/i<To2  +  CS2/1000S2  +  Cy2/co2  z  CT2/4TQ  (33) 

2 
CT2  =  DT(L)/L3  (33a) 

2 
Cs2  =  DS(L)/L3  (33b) 

_2 
Cv2  =  DV(L)/L3  (33c) 

where : 

v  =  absolute  value  of  the  medium's  mean  velocity 
across  the  acoustic  path, 

T  =  mean  absolute  temperature  (°C), 

S  =  mean  salinity  in  parts  per  thousand, 

C   =  mean  speed  sound  (m/sec), 

DT(L)  =  value  of  the  temperature  structure  function 
at  range  L, 

D  (L)  =  value  of  the  salinity  structure  function 
at  range  L, 

D  (L)  =  value  of  the  turbulence  structure  function 
at  range  L. 

The  values  for  C   which  Black  observed  were: 

-h    -1 

In  thermocline:       5.2  x  10   m 

Near  surface:         1.0  x  10  m 

-5  -1 
Below  thermocline:     2.1  x  10  ^m 
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4.   Debye  Approximation 

The  Debye  approximation  to  the  wave  equation  is  the 
most  general  model  currently  available  for  the  description 
of  the  turbulent  ocean  effects  on  sound  transmission  [Neubert 
1970].   It  contains  Rytov's  method  as  part  of  its  solution 
and  it  reduces  to  the  Born  model  when  the  single  scattering 
condition 

X   >>  2  ir  a/aU  (34) 

is  applied.   The  resultant  equation  from  the  model, 

V2  z   2   a2  k02  Le  L  (35) 

where  L  =  the  Eulerian  integral  scale  (e.g.  the  largest 
turbulent  scale),  has  no  special  long  range  nor  high  frequency 
restrictions.   In  other  words,  it  is  good  when  the  cumulative 
effects  of  small,  continuous  fluctuations  of  the  local 
turbulent  velocity  strongly  affect  phase  at  longer  ranges. 
No  use  has  been  made  of  this  model  in  the  ocean  to  date. 

C.   MODEL  APPLICATION 

Stewart  and  Grant  [1962]  conducted  a  turbulence  experiment 
near  the  surface  in  the  presence  of  large  waves.   They  showed 
the  Taylor  microscale,  a,  was  on  the  order  of  6  cm  which 
corresponds  to  the  results  (6.4  cm)  Neubert  and  Lumley  [1970] 
calculated  from  Sheehy's  [1950]  data.   Since  these  measure- 
ments correspond  well  and  fit  the  conditions  of  the  current 
NPS  experiment,  the  Neubert  and  Lumley  values  were  used  to 
provide  a  rough  basis  for  further  analysis.   These  values, 
coupled  with  experimental  parameters,  are  shown  in  Table  1. 
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The  limits  of  applicability  of  each  of  the  wave  models  and 

2 
the  associated  value  of  V  calculated  on  the  basis  of  this 

data  is  shown  in  Table  2.   It  would  appear  from  the  latter 

table  that  the  Debye  approximation  and  perhaps  the  Born 

model  for  the  wave  region  (if  the  limit  imposed  by  relation 

(6)  is  not  too  strong)  would  be  most  applicable  in  this 

2 
case.   The  large  difference  in  V  resultant  from  the  Born 

model  compared  to  that  from  the  Debye  model  could  be  the 

result  of  poor  definition  of  the  micro  and  integral  scales. 

One  should  also  note  that  the  acoustic  near-field 

criteria  for  a  piston  source  imposes  an  additional  requirement 

on  the  wavelength: 

X  >  r2/L  (36) 

where  r  =  radius  of  the  transmitter  active  disc.   Since 
relation  (36)  requires  X   >   0.6  cm,  geometrical  optics  and 
Rytov's  method  can  be  immediately  discarded. 
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Parameter 


u/C 


TABLE  1 
PARAMETER  VALUES  FOR  MODEL  TESTING 

Assumed  Value 


g 


R, 


2 

a  a 


L 
f 
A 

Cn 

u 

ng 

a 

lo 

a 

ko 

D 

r 

NL  =  Neubert  and  Lumley  [1970] 

S  =  Sheehy  [1950] 

B  =  Black  [1965] 

C  =  Calculated  from  above  data, 


1505  m/sec 

7  x  10"5 

40m 

6.4  cm 

6300 

7 

5.0  x  10"7 

200  cm 

20  khz 

7.5  cm 

-4  -1 
1.0  x  10  ni 

10.5  cm/sec 

258  cm 

2.4  cm 

2.2  cm 

4.6  x  10"4 

0.84  cm"1 

130 

10.88  cm 

Source 
Experiment 

NL 

NL 

NL 

NL 

NL 
S 
Experiment 
Experiment 
Experiment 
B 
C 
Equation  22 
Equation  21 
Equation  24 

C 
Equation  18 
Equation  27 
Experiment 


26 


TABLE  2 
RESULTS  FROM  MODEL  TESTING 


Model 


V2 


Most  Restrictive  Limit 


Geometric 


5.8  x  10  2  (19) 


X  «  1.1  x  10  2cm  (7) 


Bom  6.3  x  10  5  (13) 

(Gaussian  form) 


1.0  x  10~2«X«  15cm(12,l6) 


Bom  7.1  x  10"5  (14) 

(Exponential  form) 


1.0  x  10  2«X«  15cm(12,l6) 


Bom 

(Neubert  Limits) 


X  <  20.9cm  (18) 


Bom  (K-O)  7.1  x  10  5  (26) 


X  »  2.9  x  10  2cm  (28) 


Rytov 


1.1  x  10"2  (3D 


X  «  0.5  cm  (30) 


Debye 

t,    =  40m) 
e 


2.4  xlO1  (35) 


None 


Debye 
(Le=  3  m) 


1.8  x  10  2  (35) 


None 


HOMOGENEOUS  ASSUMPTION 


X  <<2.2cm  (6) 


ACOUSTIC  FAR-FIELD  CRITERION 


X  >  0.6  cm  (36) 


Note:  Numbers  in  parenthesis  refer  to  applicable  equation, 
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III.   EXPERIMENT  AND  DATA  ANALYSIS 

A.   EXPERIMENT 

Measurements  were  taken  from  May  2H   to  26,  1974  from 
the  Naval  Undersea  Research  and  Development  Center's  (NUC) 
oceanographic  research  tower  located  1.6  km  off  Mission 
Beach  San  Diego,  California.   The  tower  is  a  stable,  three- 
tiered  concrete  and  steel  structure  built  into  sandy  bottom 
in  water  of  18.5  m  mean  depth.   It  provides  enclosed  spaces 
for  electronic  equipment  and  open-deck  working  areas  for 
handling  instruments.   Winch  controlled  trolleys  are  avail- 
able on  three  sides  of  the  structure  on  which  devices  may 
be  lowered  to  within  1  m  of  the  bottom.   Instruments  for 
this  experiment  were  installed  on  an  aluminum  and  steel 
frame  attached  to  a  trolley  on  the  west  side  of  the  tower 
which  faced  the  predominant  swell.   A  schematic  view  of 
the  experiment  layout  is  shown  in  Figure  1. 

Basic  meteorological  and  oceanographic  data  were  avail- 
able from  equipment  permanently  installed  on  the  tower.   A 
bathythermograph  trace  obtainable  from  the  north  side  of 
the  structure  facilitated  placement  of  the  instrument  frame 
in  the  desired  vertical  position  with  respect  to  the 
thermocline . 

The  general  environmental  conditions  encountered  were 
moderate  winds ,  sea  and  swell  from  the  West.  Light  winds 
to  5  kts,  overcast  skies  and  swell  to  1  m  were  usual  in  the 
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SCHEMATIC      DIAGRAM    OF  EXPERIMENT   LAYOUT. 
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morning  (see  Table  3).   In  the  afternoon,  winds  Increased 
to  15  kts  and  seas  of  0.75  m  with  accompanying  swell  of  1  m 
was  common.   The  diurnal  variation  in  thermocline  depth  was 
approximately  12  m. 
1.   Instruments 

Temperature  measurements  were  made  with  Penwell 
Electronics  Inc.  type  k*196,  double-bead  thermistors  which 
exhibit  low  self-heating  and  well-matched,  high-resolution 
characteristics.   Calibration  of  the  thermistors,  together 
with  associated  circuitry,  was  conducted  in  the  laboratory 
under  closely  controlled,  stable  conditions.   A  Hewlett- 
Packard  HP-2  801A  quartz  thermometer  with  an  accuracy  of 
±  .006°C  was  used  as  a  temperature  standard  and  a  decade 
box  was  used  as  a  resistance  standard.   The  same  decade  box 
was  used  to  calibrate  the  thermistor  array  in  the  field  to 
eliminate  resolution  difficulties  resulting  from  the  poor 
overall  accuracy  of  the  decade  box. 

Water  particle  velocities  were  measured  using  Marsh- 
McBirney  electromagnetic  flowmeters  (models  721  and  722). 
Since  these  flowmeters  were  calibrated  by  the  manufacturer 
only  for  steady  flow  conditions,  each  probe  was  recalibrated 
for  unsteady  flow  by  oscillating  it  in  a  water  tank  at  rates 
corresponding  to  expected  experiment  velocities. 

A  Baylor  Co.  wave  staff  system,  model  13528R,  was 
used  to  measure  sea  surface  elevation.   Calibration  was 
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accomplished  at  the  NUC  tower  by  successively  shorting  out 
sections  of  the  sensor  active  element  while  measuring  output 
voltage. 

A  Naval  Research  Laboratory,  Underwater  Sound 
Reference  Division  (USRD)  type  F-33  general  purpose  uni- 
directional transducer  was  used  as  a  CW  sound  transmitter. 
It  features  two  piezoelectric  ceramic  arrays  mounted  coaxially 
which,  when  excited  in  parallel,  exhibit  exceptional  direc- 
tivity.  This  was  desired  to  reduce  interference  from  surface, 
bottom  and  equipment  frame  reflections.   The  beam  widths 
(-3db  down)  for  20kHz  and  40kHz  were  20°  and  10°  respectively. 
USRD  type  F-50  lead  zirconate-titanate  omnidirectional 
transducers  were  used  as  receivers. 

The  sonar  array  was  calibrated  in  an  anechoic  test 
tank  filled  with  fresh  water  and  the  hydrophone  calibrations 
compared  satisfactorily  with  calibrations  previously  per- 
formed at  USRD.  .  The  frequency  standard  used  was  a  Hewlett- 
Packard  HP-2801A  quartz  thermometer  configured  as  a  frequency 
meter  and  the  voltage  standards  were  an  RMS  calibrated  Fluke 
digital  voltmeter  and  a  Hewlett-Packard  model  HP  3^00A  RMS 
voltmeter.   Phase  measurements  were  conducted  using  a  Dranetz 
model  305  phase  meter.   Calibration  and  accuracy  measurements 
for  the  peripheral  equipment  associated  with  the  sound  arrays 
were  taken  concurrently  with  transducer  calibration. 

The  instrumentation  accuracy,  response  and  limitations 
are  summarized  in  Table  4. 
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2.   Instrument  Arrangement 

Two  arrays  of  thermistors  were  mounted  on  one  side 
of  the  equipment  frame  in  a  common  vertical  plane.   Seven 
thermistors  were  mounted  along  a  horizontal  member  and  7 
were  mounted  along  a  vertical  member  as  shown  in  Figure  2. 
The  instrumentation  and  abbreviated  notation  used  in  the 
text  is  given  in  Table  5.   One  thermistor  in  the  vertical 
array  was  common  to  the  horizontal  array.   The  thermistor 
probes  were  spaced  in  an  approximate  geometric  progression 
so  that  a  maximum  number  of  spatial  combinations  could  be 
achieved  for  adequate  definition  of  the  spatial  correlation 
and  structure  functions  with  a  limited  number  of  probes. 
The  minimum  spacing  (5.0  cm)  was  determined  by  the  physical 
size  of  the  thermistor  mounts.   This  was  compatible  with 
the  smallest  wavelengths  which  could  be  measured  within  the 
time  constant  of  the  thermistor.   Maximum  thermistor  spacing 
was  determined  by  physical  constraints  on  the  length  of  the 
support  pipe  to  prevent  strumming. 

Three  velocimeters  were  placed  so  that  the  horizontal 
and  vertical  particle  velocities  could  be  measured  in  the 
same  plane  as  the  thermistor  arrays.   Additionally,  they 
were  placed  in  a  vertical  line  directly  under  the  Baylor 
wave  gage  with  the  horizontal  component  measured  from  the 
west.   The  temperature  and  particle  velocity  devices  were 
offset  29  cm  from  the  support  frame  by  mounting  brackets  to 
minimize  false  turbulence  effects  generated  by  the  frame. 
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TABLE  5 
CHANNEL  ABBREVIATIONS 


Data  Channel 


Abbreviation 


Vertical  Temperatures      TMV1-TMV7 


Horizontal  Temperatures    TMH1-TMH7 


Sound  Amplitude 


Sound  Phase 


Sea  Surface  Elevation 


SAM1-SAM3 
SPH1-SPH2 
WAVE 


Horizontal  Flow  Velocity   FVAY-FVCY 


Vertical  Flow  Velocity     FVAX-FVCX 


Order  in  Fig.  2 

From  Bottom 

From  West 

From  East 

or 
From  Top 


Conductivity 


SAL1-SAL2 


From  Bottom 
From  Bottom 
From  West 
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The  sonar  array  was  separately  mounted  on  a  tripod 
base  which  pivoted  vertically  around  the  main  equipment 
frame  but  which  was  held  rigid  by  clamps  during  recording 
periods.   This  arrangement  facilitated  horizontal  and 
vertical  sound  path  measurement  (Figure  2).   The  centerline 
of  the  transmitted  sound  beam  was  57  cm  away  from  the 
vertical  plane  of  the  thermistor  arrays  to  minimize  unwanted 
reflections  from  the  other  instruments. 

Three  F-50  receivers  were  placed  along  the  center 
of  the  transmitted  beam  outside  the  near  field  at  2.0,  2.3 
and  2.9  m  from  the  transmitter  face.   The  minimum  separation 
(30  cm)  was  selected  to  maximize  the  number  of  wavelengths 
through  which  the  lowest  frequency  (20kHz)  sound  would 
travel  between  receivers  and  to  remain  near  the  lower  spatial 
scales  of  the  measured  temperature.   The  maximum  sound  path 
length  of  2.9  m  was  selected  to  minimize  reflections  from 
the  frame  and  mountings. 

3.   Pre-recording  Signal  Processing 

The  signals  from  the  wave  height  and  particle  velocity 
devices  were  recorded  directly.   All  other  signals  were 
processed  through  a  combination  of  filters,  amplifiers  and 
dc  bias  removal  devices  to  maximize  signal-to-noise  (S/N) 
ratio  before  recording. 

a.   Temperature  Circuitry 

Each  of  the  horizontal  and  vertical  temperature 
arrays  was  connected  to  a  separate  temperature  array  drawer 
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which  contained  Wheatstone  bridge  networks,  amplifiers  and 
calibration  facilities  for  seven  channels  (see  Figure  3). 
The  Wheatstone  bridge  network  provided  the  capability  to 
match  the  response  of  each  channel.   Additionally,  balance 
adjustments  were  provided  to  adjust  the  operating  point  of 
the  amplifiers  to  the  mean  temperature  so  that  distortion 
during  recording  was  minimized. 
b.   Sonar  Circuitry 

As  shown  in  the  block  diagram  (Figure  4),  signals 
from  the  F-50  transducers  were  amplified  and  then  band-pass 
filtered  at  the  operating  frequency.   The  signal  was  then 
processed  through  a  second-detector-type  demodulator.   The 
dc  bias  was  removed  from  the  detected  signal  and  the  ampli- 
tude modulation  (AM)  signal  was  then  amplified  and  low-pass 
filtered. 

Dranetz  type  305  phase  meters  were  connected 
between  the  band-pass  filter  outputs  from  receiver  hydrophones 
one  and  two  and,  similarly,  between  hydrophones  two  and  three. 
The  dc  bias  was  removed  from  the  output  of  each  phase  meter 
and  the  signal  was  then  amplified  and  low-pass  filtered. 

Before  each  run  the  dc  bias  setting  was  made  and 
the  rms  voltage  from  each  band-pass  filter  output  and  mean 
phase  shift  read  from  each  phase  meter  was  recorded.   This 
technique  allowed  much  higher  signal  amplification  to  be 
used  prior  to  recording  because  the  signal  amplifiers  were 
not  saturated  by  the  large  dc  segment  of  the  signal. 
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OCEAN      ENVIRONMENT 
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FIGURE  3. 


BASIC    TEMPERATURE    MEASUREMENT 


CIRCUIT      BLOCK    DIAGRAM.    ONLY    I  OF  THE  7  CHANNELS   IS 
SHOWN     IN  THE    TEMPERATURE    ARRAY   DRAWER     FOR   CLARITY. 
TWO     IDENTICAL    ARRAY    DRAWERS  WERE  USED    DURING    THE 
EXPERIMENT. 
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SONAR    ARRAY     IN    OCEAN  ENVIRONMENT 
2  m  
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13              1 

FIGURE  4.  BASIC    SONAR   PHASE  AND  AMPLITUDE 

MEASUREMENT     BLOCK     DIAGRAM.    TWO    PHASE    CHANNELS  AND 
THREE    AMPLITUDE     CHANNELS    WERE    USED,  HOWEVER   ONLY 
ONE     OF    EACH   IS     SHOWN     FOR    CLARITY. 
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4.   Experimental  Procedure 

Twelve  data  recording  runs  of  29  to  73  minutes 
duration  were  taken  at  various  depths  above,  in  and  below 
the  thermocline.   A  total  of  28  channels  of  data  were 
recorded  for  each  run  including  2  channels  dedicated  to 
experimental  conductivity  devices  [Kane  197*0.   Surface 
elevation,  vertical  and  horizontal  particle  velocities, 
temperature,  sound  amplitude  modulation,  sound  phase  modu- 
lation and  conductivity  were  measured  simultaneously. 

Run  depths  were  selected  on  the  basis  of  a  BT  trace 
taken  just  before  the  run,  but  the  final  run  depth  was 
chosen  to  optimize  the  degree  of  fluctuation  in  temperatures 
monitored  from  the  vertical  thermistor  array.   The  sonar 
array  was  held  in  a  fixed  position  during  each  recording 
run.   Runs  were  made  with  the  sonar  array  oriented  hori- 
zontally, vertically  and  at  an  angle  with  respect  to  the 
support  structure.   This  allowed  data  to  be  obtained  at 
various  attitudes  through  and  along  the  thermocline.   A 
summary  of  all  twelve  runs  is  given  in  Table  6. 

A  Sabre  III  and  a  Sangamo  model  3500,  both  Ik   channel 
FM  tape  recorders,  were  used  to  record  the  analog  data.   A 
timing  pulse  was  recorded  on  one  channel  of  each  recorder 
at  the  beginning  of  each  run  to  provide  temporal  continuity 
between  recorders. 

The  analog  recordings  were  digitized  using  a  hybrid 
system  including  a  Scientific  Data  System  model  5000  analog 
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computer  and  a  Xerox  Data  Systems  model  9300  digital  computer, 
Digital  data  preparation  and  analysis  was  performed  on  an 
IBM  model  360  computer.   The  digitizing  data  rate  used  was 
25Hz,  however  60Kz  electronic  noise  was  coupled  into  the 
signal,  so  low-pass  filters  were  used  to  remove  it.   Addi- 
tionally, noise  spikes  which  existed  in  some  channels  were 
removed  by  a  digital  editing  routine.   Data  handling  capacity 
was  improved  by  reducing  the  number  of  samples  by  five  using 
an  inverse  transform,  numerical  low-pass  filter  [Davidson 
1970].   The  latter  reduced  the  analysis  Nyquist  frequency 
to  2.5Hz.   The  linear  trend  and  remaining  dc  bias  were  also 
removed  before  analysis. 

Spectral  analysis  was  performed  on  26.6-  minute 
segments  of  certain  data  sets.   This  choice  dictated  a 
minimum  detectable  frequency  of  0.00625Hz.   A  standard  lag 
time  of  5$  of  the  record  length  was  chosen  which  resulted  in 
40  degrees  of  freedom  for  each  computed  spectrum.   The 
resultant  80%    confidence  limits  of  the  chi-squared  distri- 
bution were  between  0.73  and  1.30  of  the  measured  spectral 
estimates.   A  Parzen  window  was  applied  to  the  correlation 
function  during  the  spectra  calculations  to  minimize  the 
effect  of  finite  record  length  on  the  spectral  analysis. 
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IV.   EXPERIMENT  RESULTS 

During  this  experiment,  the  primary  factors  influencing 
the  results  were  the  sea  state,  the  thermocline  gradient  and 
the  position  of  the  instruments  with  respect  to  the  thermo- 
cline.  The  sea  state  remained  nearly  constant  throughout 
the  recording  periods  although  the  peak  of  the  swell  showed 
a  gradual  shift  to  higher  frequencies.   The  thermocline 
depth  and  gradient  continuously  fluctuated  under  the  influence 
of  internal  waves  which  are  characteristic  of  the  waters  off 
San  Diego  during  this  time  of  year.   In  order  to  resolve 
some  of  the  effects  of  the  thermolcine  variability,  analyses 
in  and  above  the  layer  were  required.   Experiment  runs  one 
and  two  were  chosen  for  this  purpose.   Additionally,  run 
12  was  examined  to  check  the  results  since:   1)  it  included 
data  in  and  above  the  thermocline  in  a  single  record;  2)  the 
thermocline  gradient  was  intermediate  between  the  gradient 
encountered  in  runs  one  and  two;  3)  the  sonar  array  was 
angled  toward  the  vertical  rather  than  horizontal  during  the 
run  and  4)  the  run  was  conducted  two  days  after  runs  one 
and  two  under  slightly  different  sea  conditions. 

A.   ENVIRONMENT 

The  thermocline  migrated  diurnally  from  near  the  surface 
in  the  morning  to  about  14m  in  the  evening.   It  was  generally 
diffuse  near  the  surface  but  near  the  bottom  it  was  compressed 
so  that  a  very  sharp  thermal  gradient  (2°C/m)  was  created. 
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Superimposed  upon  this  effect  was  an  apparent  four-hourly 
dilation  and  contraction  of  the  isotherms. 

Figures  5-9  show  representative  bathythermograph  (BT) 
traces  taken  periodically  during  the  experiment.   The  traces 
are  plotted  so  that  the  12°C  mark  corresponds  to  the  time 
the  BT  was  taken.   Recording  run  envelopes  in  depth  and  time 
are  included  to  show  the  position  of  the  instrument  array 
relative  to  the  thermocline.   The  11°,  12°  and  13°C  isotherms 
are  shown  to  indicate  the  relative  intensity  of  the  tempera- 
ture gradient  during  the  runs.   Figure  9  is  an  expanded 
time  plot  for  BT's  taken  during  run  12  which  shows  a  large 
internal  wave  having  a  36  min  period;  this  was  observed 
throughout  the  experiment. 

The  four-hourly  variation  in  the  thermocline  gradient 
probably  was  caused  by  internal  wave  interference  effects 
since  analysis  showed  that  these  variations  were  associated 
with  periods  of  high  temperature  variance  but  not  with  sur- 
face waves  nor  heat  budget  considerations. 

Surface  slicks  presumably  associated  with  internal  waves 
were  seen  propagating  from  the  West  at  certain  intervals 
during  the  experiment.   A  typical  example  of  the  effect  of 
these  waves  is  shown  in  the  vertical  temperature  analog 
records  of  run  number  12  (Fig.  10).   The  large  36  min  wave 
peak  passing  the  vertical  temperature  array  is  shown  by  the 
dramatic  temperature  decrease  which  occurred  at  about  1800. 
Nine  and  1.5  min  period,  small-amplitude  internal  waves 
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effects  are  seen  superimposed  upon  the  larger  internal  wave. 
The  "grass"  is  surface  swell  related.   Spectral  analysis 
showed  that  a  broad  band  of  internal  wave  periods  from 
160  sees  (maximum  detectable  by  the  method)  to  as  low  as 
30  sees  were  superimposed  on  the  large  wave  as  well. 

Figure  11  shows  the  resultant  sound  effects  from  the 
internal  wave  temperature  perturbations  of  Pig.  10  during 
run  12.   The  sound  array  was  angled  2  5°  to  the  vertical  in 
run  12.   The  sound  amplitude  modulated  signals  SAM1  and  SAM2 
were  initially  above  the  layer  and  were  not  affected  by 
the  turbulent  thermocline  losses  as  severely  as  SAM3  (below 
the  layer).   As  the  layer  moved  up,  the  interference  and 
diffraction  effects  of  the  turbulent  blobs  [Whitmarsh  and 
others  1957]  in  the  vicinity  of  SAM3  were  diminished.   This 
is  shown  by  the  increase  in  sound  pressure  level  (about 
10~3  ub)  on  the  SAM3  analog  record. 

When  the  layer  moved  into  the  vicinity  of  SAM1  and  SAM2 , 
the  sound  pressure  level  at  these  transducers  was  reduced 
by  about  6  x  10   yb  and  3  x  10   yb,  respectively.   The 
thickness  of  the  "grass"  on  these  plots  is  an  indication  of 

the  wave-induced  amplitude  scintillation  effects  (-2   to  5 

-4 
x  10   yb)  caused  by  pumping  the  turbulent  blobs  back  and 

forth  across  the  sound  path.   (Similar  results  occurred  when 

the  sound  path  was  horizontal,  but  the  scintillations  were 

then  caused  by  waves  pumping  the  thermocline  up  and  down.) 

Additionally,  the  movement  of  the  layer  (AT  a  . 3°C)  between 

SAM1  and  SAM2  caused  sound  phase  modulation,  SPH1,  to  phase 

shift  about  30°. 
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B.   TEMPORAL  STOCHASTIC  ANALYSIS 

Time  series  analyses  were  performed  on  26.67  minute 
segments  of  each  of  the  sound  channels,  the  wave  record, 
three  temperature  records  and  two  water  particle  velocity 
records.   Analyses  included  for  each  signal  the  calculation 
of  the  first  and  second  moments,  variance,  mean  period 
(calculated  from  total  zero  upcrossings) ,  autocorrelation 
function  and  autospectrum.   Also,  the  crosscorrelation, 
quadrature  and  co-spectra,  phase  difference,  coherence- 
squared  and  linear  regression  analysis  were  calculated 
between  specific  pairs  of  channels  (see  Table  7).   Examples 
of  cross-spectral  analysis  results  for  runs  two  and  12  are 
shown  in  Figs.  12-15.   Additional  results  are  given  in 
Appendix  A.   Thermistor  2  on  the  horizontal  array,  TMH2 , 
was  selected  for  temperature  comparisons  with  the  sound 
records  because  this  thermistor  was  most  central  to  the 
sound  propagation  path  during  the  runs. 

1.   Basic  Statistics 

Representative  basic  statistics  of  temperature, 
waves,  sound  amplitude  and  sound  phase  are  shown  in  Table 
8.   By  examining  the  order  of  magnitude  of  these  statistics 
in  relation  to  the  depth  of  the  run  and  the  proximity  to 
the  thermocline, certain  qualitative  relationships  between 
the  parameters  can  be  seen: 

1)   The  magnitude  of  the  temperature  variance  was 
determined  primarily  by  proximity  to  the  layer  and  its  mean 
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TABLE  7 
CHANNELS  COMPARED  IN  TEMPORAL  STOCHASTIC  ANALYSES 

Abscissas  For  Regression  Analyses 
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FIGURE    12.      Run   2   TOH2  vs   WAVE  Spectra 
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period  was  controlled  by  surface  wave  swell  until  the  layer 
was  very  close  to  the  thermocline. 

2)  Sound  amplitude  variance  was  controlled  by  proximity 
to  the  thermocline  (temperature  gradient)  when  the  sonar 
array  was  horizontal.   The  movement  of  the  thermocline 

with  respect  to  the  sound  path  must  therefore  be  the  primary 
influence  on  sound  amplitudes  in  this  configuration. 

3)  The  amplitude  variance  when  the  array  was  nearly 
vertical  was  much  smaller  than  when  it  was  horizontal. 
This  indicated  the  effect  of  turbulence  proximity  to  the 
hydrophone  was  smaller  than  the  effect  of  temperature 
gradient  motion  over  the  sound  path. 

M)   Since  the  mean  period  of  the  sound  fluctuations  was 
much  smaller  than  the  mean  period  of  the  swell  at  depth, 
turbulent  interactions  related  to  the  higher  frequency 
internal  wave  motion  of  the  thermocline  were  emphasized 
in  the  sound  measurements. 

5)  When  the  sound  array  was  vertically  oriented  vice 
horizontal,  much  larger  phase  variations  occurred  even 
though  amplitude  variations  were  smaller.   This  resulted 
from  the  large  temperature  variations  over  thepath  length 
caused  by  waves  pumping  the  thermal  gradient  up  and  down. 
This  also  showed  that  the  temperature  variation  effects 
on  sound  amplitude  and  phase  were  in  quadrature. 

6)  The  order  of  the  mean  sound  phase  variances  indicates 
that  phase  was  primarily  influenced  by  surface  waves  when  the 
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array  was  horizontal.   Since  spectral  analysis  showed  that 
the  wave  coherence  with  temperature  was  very  low  at  the 
surface  wave  frequencies,  this  effect  may  be  bubble  related. 
2.   Spectral  Analysis 

Cross-spectral  analysis  of  the  temperature,  waves 
and  sound  records  for  the  horizontal  sonar  runs  one  and 
two  showed  that  the  sound  amplitude  changes  (Figs.  13  and 
1*1)  were  in  phase  with  the  waves  and  180°  out  of  phase  with 
temperature  in  the  regions  of  highest  coherence.   This  was 
expected  in  the  stable  density  structure  encountered. 
Additionally j  the  coherence  squared  function  (simply  re- 
ferred to  as  coherence)  showed  that  only  the  long-period, 
prominent  swell  was  important  to  the  sound  amplitude  modu- 
lation because  the  higher  frequency  wave  particle  velocities 
were  much  more  rapidly  attenuated  with  depth. 

Sound  phase  (Fig.  15)  was  found  to  vary  directly 
with  temperature  at  the  internal  wave  and  surface  wave 
frequencies  although  coherence  at  the  surface  wave  frequen- 
cies was  much  smaller  than  at  the  internal  wave  frequencies. 
This  indicated  that  the  turbulence  associated  with  the  layer 
was  more  important  than  bubbles  or  surface  wave  perturbation 
of  the  thermocline  in  influencing  sound  speed  over  the 
horizontal  sound  path.   A  comparison  of  the  coherences  of 
Figs.  12  and  13  shows  that  the  temperature  effect  is  at 
least  as  important  as  bubbles  even  at  the  surface  wave 
frequencies. 
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Similar  cross-spectral  analyses  with  the  sonar  path 
angles  25°  from  the  vertical  (run  12)  yielded  markedly 
different  results.   For  the  sound  path  entirely  above  the 
layer  (SAM1),  coherence  with  temperature  was  small  over  all 
frequencies.   Small  coherence  peaks  occurred  at  0.15Hz  and 
0.7Hz  which  indicated  that  the  vertical  motion  of  the  turbu- 
lent blobs  forced  by  surface  swell  had  some  effect  in  the 
angled  sound  path. 

Sound  phase  coherence  with  temperature  was  largest 
at  the  internal  wave  frequencies  but  smaller,  broad-banded 
coherence  at  the  surface  wave  frequencies  was  also  observed. 
This  confirmed  that  the  phase  was  being  controlled  by  the 
vertical  motion  of  the  thermocline. 

Most  of  the  sound  path  associated  with  SAM2  (2.3m) 
was  not  in  the  thermocline  during  the  majority  of  run  12. 
However  the  hydrophone  was  near  the  top  of  the  layer.   As  a 
result,  its  amplitude  was  in  phase  with  temperature  at  the 
internal  wave  frequencies  and  coherence  at  other  frequencies 
was  almost  zero.   This  verified  that  the  interference  effect 
of  the  proximity  of  the  thermocline-associated  turbulence 
to  the  receiver  was  theprimary  effect  on  SAM2  variability. 

Very  high  coherence  and  180°  phase  difference  between 
TMH2  and  SAM3  at  the  internal  wave  and  the  surface  wave 
frequencies  was  observed.   Since  the  acoustic  path  associated 
with  SAM3  generally  spanned  the  thermocline,  this  indicated 
that  the  variability  of  amplitude  across  the  entire 
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thermocline  was  controlled  by  refraction  and  diffraction 
similar  to  the  horizontal  case.   The  effect  of  diffraction 
was  smaller  than  for  the  horizontal  case,  because  the  len- 
ticular shape  of  the  turbulent  blobs  present  broader  profiles 
to  the  vertical  sound  path.   The  reduced  diffraction  effects 
caused  the  variance  in  sound  amplitude  for  run  12  to  be  an 
order-of-magnitude  less  than  that  for  run  1. 

3.   Regression  Analysis 

2 

Sound  amplitude  and  phase  variances,  V  (SAM)  and 

V  (SPH),  were  plotted  versus  TMH2  variances,  V2(T),  and 

least  squares  regression  analyses  were  applied  to  the  data. 

2  2 

Linear  relationships  between  V  (SAM)  and  V  (SPH)  versus 

p 

V  (T)  were  found  to  exist  for  runs  one  and  two  but  not  for 

run  12.   This  was  because  of  insufficient  variability  in 
the  portion  of  run  12  analyzed.   The  empirical  relationship 
below  was  developed  from  the  mean  regressions  calculated 
for  runs  one  and  two: 


V2(SAM)  =  (0.3^yb2/°C2)V2(T)e"L/x  (37) 


where  L  =  acoustic  path  length  and  x  =  the  distance  required 
to  reduce  V  (SAM)  by  1/e  (determined  empirically).   x  is 
near  the  size  of  the  Taylor  microscale. 

The  derivation  of  Equation  (37),  plots  of  the 
variances  over  the  run  segments  and  the  regression  plots 
are  given  in  Appendix  B. 
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C.   SPATIAL  ANALYSIS 

Temperature  spatial  correlation  and  structure  functions 
were  calculated  directly  from  ensemble  averages  of  26.67  min 
records  from  all  possible  combinations  of  thermistor  pairs 
along  the  vertical  thermistor  array.   Twenty-two  spatial 
lags  with  values  between  0  and  260  cm  were  thus  determined 
from  seven  thermistors.   The  value  of  CT  (see  Equation  33a) 
and  the  Taylor  microscale  as  well  as  the  Eulerian  integral 
scale  were  derived  from  the  above  functions  (Table  9). 
However,  since  large  variations  were  seen  in  these  functions 
the  validity  of  the  calculated  scales  was  questionable. 

In  an  attempt  to  remove  the  non-stationary  effects 
resulting  from  internal  waves,  analyses  using  ensemble 
averages  of  short  segments  having  nearly  constant  turbulent 
intensity  were  conducted  for  runs  two  and  twelve.   The 
segments  were  chosen  so  that  examples  of  low,  medium  and  high 
turbulent  intensity  could  be  examined  (Table  10).   Because 
the  thermocline  steepness  was  highly  variable  along  the 
vertical  thermistor  array,  however,  only  the  medium  inten- 
sity correlation  function  was  smooth  enough  to  provide  satis- 
factory Eulerian  integral  scales  and  the  structure  function 

of  all  samples  varied  widely.   Values  of  the  Taylor  microscale, 

2 

Eulerian  integral  scale  and  C„   generated  from  hand-smoothed 

results  are  given  in  Table  11.   The  correlation  function  and 
structure  function  values  as  well  as  the  mean  variance  for 
the  vertical  array  thermistors  are  given  in  Appendix  C. 
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TABLE    9 
Temperature   Spatial   Scales    (cm)    for  26.67  mln  Ensembles 


Run  # 

Taylor 
Micro scale 

15 

Taylor 

Microscale 

40 

Eulerian 
Integral  Scale 

1 

36 

2 

178 

158 

175 

12 

60 

100 

36 

V 


1.45x10 


-4 


3.6x10 
8.3x10 


,-3 
,-3 


TABLE   10 
Relative  Turbulent   Intensity   for   Special  Spatial  Analyses 


Qualitative 
Intensity 

Ensemble 
Size  (min) 

TMH2  Mean 
Variance (°C)2 

Range 
(°C)2 

Run 

Low 

4.67 

0.27xl0-1 

±0.77xl0"2 

12 

Medium 

4.67 

4.24X10"1 

±4.95xl0_1 

2 

High 

3.73 

7.79X10"1 

±2.23xl0~1 

2 

TABLE    11 
Temperature   Spatial  Scales    (cm)    for  Special  Ensembles 


Qualitative 
Intensity 

Taylorx 
Microscale 

23 

Taylor 
Microscale 

6.8 

Eulerian 
Integral  Scale 

24 

C  2 

LT 

Low 

3.3xlO-3 

Medium 

71 

40 

35 

4.4xlO"3 

High 

15 

20 

63 

5.4xl0"3 

Taylor    =  Taylor  microscale  from  spatial  correlation 

p 

Taylor  =  Taylor  microscale  from  structure  function 
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D.   APPLICATION  OF  MODELS 

While  it  is  apparent  that  some  of  the  scale  estimates 
calculated  from  the  vertical  thermistor  array  may  be  erron- 
eous, it  is  instructive  to  apply  the  theoretical  models 
for  the  sound  amplitude  variance  to  see  which  model  is 
most  applicable  In  the  non-stationary  case.   Table  12  lists 
the  equations  necessary  to  calculate  the  basic  parameters 
for  application  of  the  models.   For  the  purposes  of  this 
examination,  the  values  of  the  Richardson  number  (Rt)>  the 

Prandtl  number  (a,  )  and  u/c   were  assumed  to  be  the  same 

k         o 

as  reported  by  Neubert  [1970]  (see  Table  13  for  values). 
The  speed  of  sound  is  used  to  calculated  Wilson's  Equation 
[Urick  1967]. 

Tables  13  and  1^  list  the  resultant  parameters  of  each 
sample  for  which  scales  were  determined.   It  is  noteworthy 
that  the  values  of  the  refractive  index  (a)  calculated  above 
are  an  order  of  magnitude  above  the  near  surface  results 
reported  by  other  investigators  [Haley  1972,  Medwin  1973] 
(see  Table  15)  except  in  the  low  turbulent  intensity  cases. 
This  is  further  evidence  that  the  relative  influence  of  the 
turbulent  thermocline  in  sound  scattering  is  much  greater 
than  bubbles. 

If  the  order  of  magnitude  of  the  Taylor  microscale  is 

correct  (see  Table  13  for  comparative  values)  the  values  of 

2 

a  a  are  one  to  four  orders  of  magnitude  greater  than  previously 
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TABLE  12 
EQUATIONS  USED  FOR  MODEL  PARAMETER  EVALUATION 

Parameter  Equation 

cQ  cQ  =  1449  +  4.6T  -  0.55T2  +  0.17d 


a2  =  [4.62V2(T)  +  0.0552v\t)]/co2 


Xg  Ag  =  a/ak  =  2.65a 


ng       .        ng  -  Ag4(Rt/i5)  =  2.07  x  10V 


Cn  Cn2  "  S'/^o2 


lQ*  l0*   =   0(2. 2/2. 4a)  =  0(0. 9a) 


T     =  mean  temperature  (°C) 

d     =  depth  (m) 

2  2 

V  (T)  =  variance  in  temperature  (°C  ) 

a     =  Taylor  microscale  (cm) 

n     =  Kolmogorov  microscale  (cm) 


*  equation  for  I      assumes  u/c   is  a  constant 

o  o 
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TABLE  13 
Experiment  Parameter  Values  for  Model  Tests 


Parameter 

RUN  1 

RUN  2 

RUN  12 

cQ(m/sec) 

1516 

1498 

1515 

a 

4.5xl0"4 

2.1xl0"3 

-4 
7.9x10 

L  (cm) 

200 

200 

200 

a  (cm) 

15 

178 

60 

c 

13.5 

160 

5* 

X    (cm) 

7.5 

7.5 

7.5 

kQ  (cm-1) 

0.84 

0.84 

0.84 

Cn 

-4 
1.45x10 

1.05x10"^ 

-4 
2.4x10  4 

Le  (cm) 

36 

175 

36 

ak 

7 

7 

7 

X      (cm) 

40 

472 

159 

ng  (m) 

10.5 

21.4 

19.2 

Rt 

6300 

6300 

6300 

2 
a  a 

3.0xl0"6 

-4 
7.9x10 

3.7xl0"5 
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TABLE  14 
Short  Sample  Parameter  Values  for  Model  Tests 


Parameter 

Low 

Intensity 

Medium 
Intensity 

High 
Intensity 

c   (m/sec) 

1512 

1508 

1500 

a 

5.0x10' 

-4 

2.0x10" 

•3 

2.7xl0~3 

L  (cm) 

200 

200 

2  00 

a  (cm) 

23 

71 

15 

* 

21 

64 

13 

X  (cm) 

7.5 

7.5 

7.5 

kQ  (cm-1) 

0.84 

0.84 

0.84 

Cn 

1.05x10"^ 

l.l6xlC 

r* 

1.29xl0"4 

Le  (cm) 

24 

35 

63 

ak 

7 

7 

7 

Xg  (cm) 

61 

188 

40 

n   (cm) 

2.8 

8.5 

1.8 

Rt 

6300 

6300 

6300 

a  a  (cm) 

5.8x10" 

-6 

2.8x10" 

-4 

l.lxl0-Z| 
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reported.   This  result  has  a  direct  relationship  (in  the 
Born  model)  to  the  variance  of  sound  amplitude. 

The  variances  of  sound  amplitude  calculated  by  each 
model  previously  discussed  (Geometric,  Born,  Rytov,  Debye) 
are  presented  in  Tables  16  -  19.   In  general,  the  best  fit 
to  the  measured  sound  variances  were  calculated  from  the 
Born  model  using  the  Krasil 'nikov-Obhukov  formulation.   A 
surprisingly  good  fit  to  the  true  variance  was  also  obtained 
from  the  geometric  model  in  the  mid-intensity  range.   Because 
of  the  above,  the  stimates  of  the  Taylor  microscale  seem 

reasonable. 

2 
Since  the  accuracy  of  the  scales,  C   and  L  ,  the  Eulerian 

integral  scale,  is  in  question,  no  firm  statement  can  be 
made  about  the  adequacy  of  the  Rytov  and  Debye  models. 
However,  because  L  is  thought  to  be  correct  for  the  mid- 
intensity  ensemble,  the  Debye  model  appears  incorrect. 
Presumably,  calculation  of  the  structure  function  and  corre- 
lation function  from  the  thermistors  in  the  horizontal  array 
over  stationary  ensembles  would  resolve  this  question. 
Unfortunately,  this  was  not  accomplished  in  time  for 
inclusion  in  this  thesis. 
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TABLE  16 
Geometric  Model  Experimental  Results 
V2  =  i|/15  /F  a2(L/a)3     ;    7.5  «  £  2/L   (cm) 


Ensemble 

Model  V2 

Observed  V2 

(ub)2 

£02/L   (cm) 

Run  1 

4.8xlO_Z| 

1.12xl0~4 

0.91 

Run  2 

3xl0"6 

1.45xl0"3 

128 

Run  12 

l.lxlO"5 

4.05xl0~5 

146 

Low  Intens 

ity 

7.8xl0"5 

7.1xl0"6 

2.2 

Mid  Intens 

ilty 

4.2xl0~5 

9.5xl0"5 

, 

21 

High  Intensity 

8.2xl0"3 

-4 
3.0x10 

0.85 

TABLE  17 

Born  Model  Experimental  Results 
[Krasil'nikov  and  Obhukov] 

V2  =  k02a2aL(l  -  1/Dtan"1  D)   ;     D  =  TTL/kQa 

?  2 

Ensemble  Model  V     Observed  V  D 


Run  1  4.3xl0_Ij  1.12xl0-i|  50 

Run  2  7.9xlO_2<  1.45xl0~3  4.2 

Run  12  3.74xl0~5  4.05xl0~5  12.5 

Low  Intensity  5.8xl0~6  7.1xl0-6  32.5 

Mid  Intensity  2.8xl0-i|  9.5xl0~5  10.5 

High  Intensity  1.1x10"^  3.0xl0_i|  50 

D  >>  1   -*■  Wave  Region 
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TABLE  18 
Rytov  Model  Experimental  Results 


V2  =  161  Cn2  L  11/6  -7/6 


Ensemble       Model  V2 


Run  1  5.2xl0"3 

Run  2  2.8xl0~3 

Run  12  1.5xl0"2 

Low  Intensity  2.8xl0**3 

Mid  Intensity  3.4xl0~3 

High  Intensity  4.2xl0~3 


7.5   <<  3/l0i 

/L 

•a 

Observed  V 

/l0   /L 

1.12xl0~4 

5.5 

1.45xl0~3 

150 

4.05xl0~5 

35 

7.1xl0"6 

9.9 

9.5xlO"5 

44 

-4 

3.0x10 

5.2 

TABLE  19 
Debye  Model  Experimental  Results 


2              2      2 

V     =   2aVLeL 

Ensemble 

Model  V2 

2 
Observed  V 

Run   1 

2.1xl0~3 

1.12X10"*1 

Run   2 

0.22 

1.45xl0~3 

Run   12 

6. 3xl0~3 

4.05xl0~5 

Low   Intensity 

1.7xlO"3 

7.1xl0~6 

Mid   Intensity 

4xl0"~2 

9.5xl0~5 

High  Intensity 

0.13 

3.0xl0_!| 
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V.   CONCLUSIONS 

The  theories  of  turbulence  which  have  been  developed 
for  the  ocean  require  some  degree  of  stationarity  in  the 
medium  before  they  may  be  applied.   This  paper  has  examined 
four  basic  stochastic  models  of  acoustic  wave  propagation 
under  the  most  severe  conditions  of  non-stationarity ,  e.g. 
near  a  sharp  thermocline  which  is  being  oscillated  by  both 
large  surface  swell  and  by  internal  waves. 

The  experimental  conditions  specified  that  the  solution 
form  would  be  in  the  wave  region  for  which  the  Born  model 
would  be  most  accurate.   The  solution  of  the  Krasil *nikov- 
Obhukov  form  of  the  Born  model  from  directly  measured, 
ocean  parameters  gave  reasonable  results. 

The  solution  of  the  Rytov  and  Debye  models  for  the  same 
conditions  yielded  results  two  orders  of  magnitude  higher 
than  was  actually  measured.   Because  evaluation  of  scales 
was  difficult,  this  was  expected  even  though  both  models 
are  more  general  than  the  Born  model  and,  in  fact,  contain 
its  solution  in  a  subset.   This  also  showed  the  determination 

of  the  Taylor  microscale  was  more  reliable  than  the  deter- 

2 

mination  of  the  structure  constant,  C„  ,  or  the  Eulerian 

integral  scale  for  measurements  across  the  thermocline. 

An  attempt  was  made  to  apply  the  correlation  function 
and  structure  function  piecewise  to  small  sections  of  data 
which  were  thought  to  be  relatively  stationary.   For 
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application  across  the  thermocline,  this  method  did  not  work 
well  because  selection  of  such  sections  is  very  difficult. 

The  results  of  the  time  series  analyses  showed  that  the 
losses  over  the  acoustic  path  were  greatest  when  the  array 
was  horizontal  near  the  thermocline  because  the  diffractive 
losses  were  an  order  of  magnitude  greater  than  those  for 
the  nearly  vertical  case.   Interference  effects  due  to  the 
proximity  of  turbulence  to  the  receiving  hydrophone  accounted 
for  a  0.l4fS  signal  loss  and  scintillations  due  to  surface 
swell  accounted  for  0.07%  signal  modulation.   The  effects 
of  bubbles  were  small  because  the  measurements  were  rela- 
tively deep  and  the  acoustic  frequency  was  much  lower  than 
the  peak  bubble  resonant  frequency.   Measured  loss  over  the 
2m  path  was  on  the  order  of  5%   and  phase  shifts  of  up  to 
30°  occurred  for  these  worst-case  conditions. 

Current  sound  propagation  models  based  upon  empirical 
data  have  undoubtedly  included  these  effects  in  surface 
scattering  or  reverberation.   But  the  roughness  description 
for  the  sea  surface  used  by  sonar  models  is  based  upon  sea 
state.   The  losses  here  described  are,  instead,  related  to 
the  swell  wave-induced  water  particle  velocities  and  internal 
waves. 
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APPENDIX  A 

Spectral  Analysis  Plots 

This  appendix  includes  examples  of  cross-spectral 
analysis  plots  so  that  qualitative  comparisons  between 
sound,  temperature  and  waves  during  the  experiment  can  be 
made.   Included  are  spectra  for  1)  sound  amplitude  versus 
temperature,   2)  sound  phase  versus  temperature,   3)  wave- 
related  effects  and   *0  comparison  of  one  sound  channel 
with  another  for  the  same  run.   The  plots  included  are 
listed  below. 

Sound  Amplitude  Versus  Temperature 

Run  #  Channels  Compared 

1  SAM1  -  TMH2 

1  SAM2  -  TMH2 

1  SAM3  -  TMH2 

2  SAM1  -  TMH2   (see  text) 
2  SAM2  -  TMH2 

2  SAM3  -  TMH2 

12  SAM1  -  TMH2 

12  SAM2  -  TMH2 

12  SAM3  -  TMH2 

Sound  Phase  Versus  Temperature 

Run  #  Channels  Compared 

1  SPH1  -  TMH2 

2  SPH1  -  TMH2   (see  text) 
2  SPH2  -  TMH2 

12  SPH1  -  TMH2 
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Wave-Related  Effects 


Run  # 

Channels  Compared 

1 

TMH2  -  WAVE 

1 

SPH1  -  WAVE 

•   2 

SAM1  -  WAVE 

6 

TMH2  -  FVBY 

12 

TMH2  -  WAVE 

Sound  Channel  Comparisons 

Run  #  Channels  Compared 

2  SAM1  -  SAM2 

2  SAM1  -  SAM3 

2  SAM2  -  SAM3 

2  SPH1  -  SPH2 
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APPENDIX  B 
Analyses  of  Variance  Regressions 

Variances  of  56  sec  samples  of  data  which  included  1^00 
data  points  each  were  calculated  along  26.67  min  segments 
of  certain  channels  of  runs  one,  two  and  12.   Regression 
analyses  for  sound  amplitude  and  phase  versus  temperature 
were  computed  using  these  variances.   The  resultant  plots 
of  the  variances  with  the  regression  curves  are  shown  after 
Table  B5.   A  summary  of  the  regression  equations  is  given 
in  Table  Bl.   The  following  mean  variance  equations  were 
computed  from  the  average  of  the  regression  equations  listed 
for  runs  one  and  two: 


V2(SAM1)  =  3.2  x  10"3  V2(T)  (Bl) 

V2(SAM2)  =  1.6  x  10~3  V2(T)  (B2) 

V2(SAM3)  =  4.1  x  lO"4  V2(T)  (B3) 


2 

If  V  (T)  is  eliminated  from  the  above  equations  (stationary 

assumption),  then  because  SAM1  and  SAM2  are  30  cm  apart  and 
60  cm  separates  SAM2  and  SAM3: 


V2(SAM)  cC  (|)~L/3°  cm  (B4) 


where  L  =  sound  path  length, 
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This  equation  can  be  checked  by  comparing  the  sound 
amplitude  variances  for  specific  56  sec  samples  (where  the 
sample  is  stationary)  and  similarly  deriving  the  relation- 
ship 


V2(SAM)  oc  (|)"L/X  (B5) 


2 

where  x  is  the  distance  required  to  reduce  V  (SAM)  by  0.5. 

This  is  done  in  Tables  B3  -  B5 .   The  following  equations 
result : 


P        n-L/34.5 
Run  1:   V  (SAM)oc  |        cm  (B6) 


Run  2:   V2(SAM)<^  ^~h/^2    cm  (B?) 


Conversion  from  base  ~-  to  base  e  and  placing  the  path 
length  reference  at  lm  vice  2m  yields  (from  equations  Bl, 
B4,  and  B5) : 


V2(SAM)  =  (0.3^yb2/°c2)V2(T)e"L/x        (B8) 


where : 

L  =  acoustic  path  length  in  meters, 

x  =  temperature  spatial  scale  for  which 
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V2(SAM)  =  V2(SAM)    /e 


2 

V  (T)  =  temperature  variance 


p 

V  (SAM)  =  variance  of  sound  amplitude. 


And  the  spatial  scales  in  (B4)a  (B6)  and  (B7)  become  43cm, 

50cm  and  6lcm,  respectively. 

2  2 

Since  V  (SAM)  is  known  to  be  a  function  of  a  ,  the 

variance  of  the  index  of  refraction,  and  a  =  2.1  x  10 

for  run  two  (which  had  the  best  regression  fits) : 


V2(SAM)  =  7.7  x  10'  a2  V2(T)  e~L/x       (B9) 


Assuming  that  x  <<  L,  the  following  approximation  can  be 
made: 


e-L/x  a   x  s  x  (B10) 

X+J-i     Li 


Substituting  Equation  B(10)  into  Equation  B(9)  and  multl- 

p 
plying  by  L   to  compensate  for  the  spherical  spreading  loss 

yields: 


V2(SAM)  =  7.7  x  10   a2  V2(T)  xL         (Bll) 
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which  is  similar  to  the  exponential  form  of  the  Born 
approximation: 


V2  =  k  2  a2  aL  (B12) 


if  k   is  similarly  included  in  Equation  (Bll): 


V2(SAM)  =  1.1  x  10  5  k  2  a2  V2(T)  xL.    (B13) 


The  plots  following  Table  B5  are  included  to  allow 
qualitative  comparisons  of  conditions  between  runs  and  of 
the  variances  between  channels  during  the  same  run.   For 
each  run  analyzed  the  variance  versus  time  plots  are  given 
followed  by  regression  plots  between  various  channels. 
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TABLE  Bl.   Regression  Equations  for  Sound  Versus  Temperature 
Variances 

AMPLITUDE 

Channel     Run#  Equation 

SAM1         1  V2(SAM1)  =  E3.^V2(T)]  x  10~3 

2  V2(SAM1)  =  [3.0V2(T)]  x  10"3 

12  V2(SAM1)  -  [4  -  2V2(T)]  x  lO-*1 


SAM2 


1 

2 

12 


V2(SAM2)  =  [2.0V2(T)]  x  10~3 
V2(SAM2)  =  [1.2V2(T)]  x  10~3 
V2(SAM2)  =  [0.7  +  1.5V2(T)]  x  10"5 


SAM3 


1 

2 

12 


V2(SAM3)  =  [4.8V2(T)]  x  10"21 

V2(SAM3)  =  E3.3V2(T)]  x  10"11 

V2(SAM3)  -  [0.2  +  4V2(T)]  x  10_Z< 
PHASE 


Channel 

Run# 

SPH1 

1 

2 

12 

SPH2 

1 

2 

12 

Equation 
V2(SPH1)  =  1  +  i)lV2(T) 
V2(SPH1)  =  0.2  +  0.7V2(T) 
V2(SPH1)  =  3.7  +  77V2(T) 


V2(SPH2)  =  2  +  8V2(T)* 
V2(SPH2)  =  2  +  22V2(T) 
V2(SPH2)  =  2  -  25V2(T) 


*From  20  minute  record. 
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TABLE  B2.      Run   1  Variances 


sample 
variance 


Run  1 
26  2/3  min 


Run  1-18 
56  sec 


Run  1-13 
56  sec 


Run  1-25 
56  sec 


vT2(°c2) 
vSAm(VIb2) 

VSAM2(^ 
VSAM3(^2) 


2.2  x  10 


-2 


-5 


11.0  x  10 


7.6  x  10"5 


2.5  x  10 


-5 


16.1  x  10~2  3.2  x  10"2  0.43  x  10"2 

52.0  x  10~5  12.0  x  10~5  18  x  10"5 

39.0  x  10~5  7.2  x  10"5  13  x  10~5 

9  x  10~5  3.3  x  10"5  3.2  x  10"5 


TABLE   B3.      Sound  Amplitude   Distance   Loss   Factors    (Run   1) 


sample 

Run  1* 

Run  1-1 

Run  1-2 

Run  1-3 

distance 

26  2/3  min 

56  sec 

56  sec 

56  sec 

30  cm 

1.H5 

1.33 

1.67 

1.38 

60  cm 

3.04 

4.33 

2.18 

4.06 

90  cm 

4.40 

5.78 

3.64 

5.63 

quantity 
distance 

mean 
loss 
factors 

standard 
deviation 
in  loss 
factors 

norma  li  zed 

loss 

factors 

accepted 

loss 

factors 

30  cm 

1.46 

0.18 

1.00 

1.15 

60  cm 

3.52 

1.17 

2.41 

2.  30 

90  cm 

5.02 

1.19 

3.44 

3.45 

*not    included   in  mean   loss    factors 
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TABLE   B4. 

Run   2  Variances 

Sample 

Run  2 

Run  2-9 

Run  2-26 

Run  2-20 

variance 

26  2/3  min 

56  sec 

56  sec 

56  sec 

vT2(°c2) 

46.2  x  10~2 

98.1  x  10"2 

143.4  x  10"2 

121.4  x  10"2 

Vsm(^ 

140  x  10"5 

311  x  10"5 

422  x  10"5 

420  x  10"5 

VSAM2(^2) 

56.0  x  10"5 

122  x  10"5 

146  x  10"5 

178  x  10~5 

VSAM3(yb2) 

17.0  x  10"5 

30.5  x  10~5 

40.3  x  10~5 

46JD  x  10"5 

TABLE  B5.   Sound  Amplitude  Loss  Factors  with  Distance  (Run  2) 


sample 

Run  2* 

Run  2-9 

Run  2-26 

Run  2-20 

distance 

26  2/3  min 

56  sec 

56  sec 

56  sec 

30   cm 

2.5 

2.5 

2.9 

2.4 

60    cm 

3.29 

4.0 

3.62 

3.9 

90   cm 

8.24 

10.2 

10.47 

9.1 

quantity 
distance 

mean 
loss 
factors 

standard 
deviation 
in  loss 
factors 

normalized 

loss 

factors 

accepted 

loss 

factors 

30    cm 

2.6 

0.3 

1.0 

1.4 

60   cm 

3.8 

0.20 

1.5 

2.8 

90   cm 

9.9 

0.73 

3.8 

4.2 

*not  included  in  mean  loss  factors. 
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^SCRLF-LOOE-Ql   UNITS    INCH..  =  VMtmH2)    C0^) 

URR1RNCE    FIT 
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Run!  i2~ 


K-SCRLE----2.00E+02  UNITS  INCH- 
Y-SCRLE----5.00E-06  UNITS  INCH. 

URR1RNCE   FIT 


-    T"v\e    C  Sec) 

=    \lz   CSM^I    )     (VOLTS*") 
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RurO  \X 


K-SCRLf--2.G0E>02  UH11S  INCH. 
V-SCRLE^S.00E-07  UNITS  INCH 

UnRIRNCE  FIT 


-time:  (sec/) 

a  Vx(S«Wl)  (  VOLTS*) 
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K-SCRLE-2. 00E+02   UNITS    INCH.  =   t/m£  C^ec) 
Y-SCRLE-5.00E-06  UNITS    INCH.-    vMsfiv\3)    (votrs2") 

URRIflNCE  FIT 
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K-SCRLF«2-"00E+02   UNITS    INCH.  -  t.aae  (sec) 
"-SCRLF---1.0QE+01    UNITS    INCH.-  v1  (  sp»0  C  de^ees  ) 

URRIRNCE   FIT 

132* 


OJS 


Ou1^ 


CiO 


RmsJ   IX 


K-SCRLF---2.00F.+02  UNITS    INCH.    =  time(sec) 
V-5CRLF---5.-00E+00  UNITS    INCH.   =   VM*PH*0    Ueyees  ; 

URRIRNCE  FIT 
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+ 


+ 


+ 


+■ 


002 


not 


OO'i 


RUM  IX 


K-SCRLL---1.00E-01  UNITS    INCH.  =  V*CtmHx)  (V) 
y-SCRLE-5.00E.-06  UNITS    INCH.   -    \ZM*A*U  <^6UT*  ) 

REGRESSION      PLOT 
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027 


+  + 

+ 


Q3'2 


00  5 


3C8 


037 


RUM  IX 


X-SCRLE--1. 001-01   UNITS  INCH- 
V-SCRIX--5.  001  -0?  UNITS  INCH. 

REGRESSION     PLQ 


=■  V^iS^Ml)  (VOLTS*) 
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Runj  IX 

■  ■  ■■--» 

K-5CRL T-l  OOF  -01  UNJTS  INCH.  = 
V-SCRLF.-5.00E-06  UNJTS  INCH.  = 

REGRESSION     PLOT 


V2-  (SPw,3)  (  VOUTS1) 


138 


RUK)    IX 


K-SCRLE--=1.00E-01   UNITS   INCH. 
V-SCRLE>I. OOE+01   UNITS    INCH. 

REGRESSION     PLOT 


--   V1"  (  SPBI  )   (  decree  $>*") 
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8  . 


++ 


+ 
+ 


+ 


022 


32: 


Run)  ix 


K-SCRLF.-1.00E-01  UNITS  INCH 
V-SCRLF.----5.00E-02  UNITS  INCH 


REGRESSION     PLOT 


Vx  (  SPH2)  (de^ree^) 


1  J<0 


++ 

■h  + 


+ 


+ 


+ 


ooc 


oc? 


ZZ-J 


oz: 


::~ 


Run)     IX 


K-SCRL 


F----1.00F-01   UMTS    INCH.  =  ^(wfiv£)   C  ^ 
02  UNITS   1NCH.=    V1"  (tmhs)  (»&  ) 


V-SCRLF.----5.00E 


REGRESSION     PLOT 
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APPENDIX  C 
Spatial  Correlation  and  Structure  "Functions 

The  values  of  the  spatial  correlation  function  and 
structure  function  calculated  during  spatial  stochastic 
analyses  are  tabulated  in  Tables  CI  -  C6.   The  functions 
were  calculated  along  the  vertical  thermistor  array. 
The  variance  of  the  temperature  at  each  thermistor  over 
the  ensemble  is  also  given. 

The  Tables  include: 


Ensemble 
Run  1,  26.67  min 
Run  2,  26.6"  min 
Run  3,  26.67  min 
Low  turbulence  intensity 
Middle  turbulence  intensity 
High  turbulence  intensity 


Table  # 

CI 

C2 

C3 

Cl 

C5 

C6 

112 


TABLE  CI.   Run  1,  Vertical  Spatial  Correlation  and 
Structure  Function. 


Thermistor  No. 

Variance 

(°c2) 

1 
2 

3 
i\ 

5 
6 
7 

0.5999 
1-0970 

0.0540 
0.3090 
0.2722 
0.1359 
0.0747 

Correlation 
Function 

Structure 

Function 

Lag 

1.0000 

0.0 

0.0 

0.9635 

0.2081 

5.0000 

0.3745 

0.5841 

15.0000 

0.3654 

0.4673 

20.0000 

0.5507 

0.3978 

25.0000 

2.6243 

0.9860 

25.0000 

0.5208 

0.5926 

30.0000 

0.1736 

0.3503 

45.0000 

0.2716 

0.2302 

55.0000 

0.2668 

0.7407 

80.0000 

0.2430 

0.9232 

85.0000 

0.0721 

0.3287 

100.0000 

0.5209 

3.1869 

135.0000 

0.9274 

3.3103 

150.0000 

0.9724 

3.0852 

155.0000 

0.3622 

1.6778 

160.0000 

0.6725 

1.9936 

175.0000 

0.6874 

1.8289 

180.0000 

0.5094 

3.5094 

180.0000 

0. 3640 

1.9750 

205.0000 

0.1944 

3.9140 

235.0000 

0.1514 

2.1753 

260.0000 
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TABLE  C2.   Run  2,  Vertical  Spatial  Correlation  and 
Structure  Function 


Thermistor  No. 

Variance 

(°C2) 

1 
2 

3 
4 

5 
6 

7 

0.2154 
0.1750 
0.3544 
0.5783 
0.3910 
0.4338 
0.5170 

Correlation 
Function 

Structure 
Function 

0.0 

Lag 

1.0000 

0.0000 

1.0592 

0.3607 

5.0000 

1.0090 

0.3676 

15.0000 

0.9452 

0.0160 

20.0000 

1.0354 

0.0388 

25.0000 

1.0891 

0.4103 

30.0000 

0.9706 

0.0749 

45.0000 

1.1503 

0.1313 

55.0000 

1.0532 

0.2159 

80.0000 

1.0844 

0.6326 

85.0000 

0.9813 

0.2662 

100.0000 

0.5463 

0.2614 

135.0000 

0.6164 

0.6939 

150.0000 

0.5559 

0.3357 

155.0000 

0.5010 

0.4553 

160.0000 

0.5551 

0.9197 

175.0000 

0.4985 

0.5539 

180.0000 

0.5580 

0.4413 

180.0000 

0.4878 

0.6848 

205.0000 

0.5458 

O.6783 

235.0000 

0.4557 

0.9619 

260.0000 
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TABLE  C3.   Run  12,  Vertical  Spatial  Correlation  and 
Structure  Function 


Thermistor  No. 

Variance  (°C2) 

- 

1 
2 

3 
4 

0.2792 

0.2433 
0.1690 
0.1444 

5 
6 

7 

0.1275 
0.0743 
0.0083 

Correlation 

Structure 

Lag 

Function 

Function 
0.0 

1.0000 

0.0 

1.0558 

0.0211 

5.0000 

1.1642 

0.1301 

15.0000 

1.0859 

0.1540 

20.0000 

0.7096 

0.1640 

25.0000 

1.5349 

1.0283 

25.0000 

0.7463 

0.2231 

30.0000 

0.7552 

0.3984 

45.0000 

0.1447 

0.3588 

55.0000 

0.1852 

0.6946 

80.0000 

0.1969 

0.8039 

85.0000 

0.1995 

0.9920 

100.0000 

0.8057 

1.6228 

135.0000 

0.6778 

1.6854 

150.0000 

0.6295 

1.6494 

155.0000 

0.6143 

2.2873 

160.0000 

0.5148 

2.2931 

175.0000 

0.4773 

2.2355 

180.0000 

0.4121 

I.6676 

180.0000 

0.2998 

2.1737 

205.0000 

0.1487 

1.6762 

235.0000 

0.1214 

2.0123 

260.0000 
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TABLE  C4.   Low-Intensity,  Vertical  Spatial  Correlation 
and  Structure  Function 


Thermistor  No. 

Variance  (°C2) 

1 
2 

3 
4 

5 
6 

7 

0.0666 
0.0774 
0.0062 
O.OO38 
0.0028 
0.0114 
0.0014 

Correlation 
Function 

Structure 

Function 

0.0 

Lag 

1.0000 

0.0 

1.0785 

0.2173 

5.0000 

1.3019 

O.986O 

15.0000 

1.0129 

1.1896 

20.0000 

0.2793 

4.5337 

25.0000 

9.4721 

32.6599 

25.0000 

0.2274 

5'.  0118 

30.0000 

0.0791 

6.1497 

45.0000 

0.1556 

4.2863 

55.0000 

0.4158 

0.6737 

80.0000 

0.4850 

0.9095 

85.0000 

0.6269 

1.4669 

100.0000 

3.1748 

23.5975 

135.0000 

1.6307 

25.8445 

150.0000 

1.5337 

25.6640 

155.0000 

2.9846 

20.1207 

160.0000 

1.5724 

22.1038 

175.0000 

1.1432 

22.5880 

180.0000 

0.7332 

30.3575 

180.0000 

-1.6798 

31.3264 

205.0000 

0.5584 

27.1199 

235.0000 

1.1856 

22.0083 

260.0000 
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TABLE  C5.   Mid-Intensity,  Vertical  Spatial  Correlation 
and  Structure  Function 

2 
Thermistor  No.  Variance  (°C  ) 


1 

2 

3 
4 

5 
6 

7 

0.3603  x 
0.3726  x 
0.8536  X 
0.1329 

0.1011 
0.8996  X 
0.3400  X 

10" 
10' 

10" 

10" 
10' 

-5 
-3 
■1 

-1 
-1 

Correlation 
Function 

Structure 
Function 

0.0 

Lag 

1.0000 

0.0 

1.0345 

0.2462 

5.0000 

0.9265 

0.3068 

15.0000 

0.8702 

0.1042 

20.0000 

0.8345 

0.2212 

25.0000 

0.0002 

0.0034 

25.0000 

0.87^2 

0.4570 

30.0000 

0.6962 

0.3425 

45.0000 

0.2549 

0.7168 

55.0000 

0.2131 

0.9103 

80.0000 

0.2374 

1.1769 

85.0000 

0.1709 

0.8393 

100.0000 

0.0130 

0.8224 

135.0000 

0.0205 

.1.2779 

150.0000 

0.0134 

0.9769 

155.0000 

-0.0008 

0.8464 

160.0000 

-0.0005 

1.3162 

175.0000 

-0.0008 

1.0017 

180.0000 

0.0149 

0.8641 

180.0000 

-0.0004 

0.8911 

205.0000 

0.0154 

0.3093 

235.0000 

0.0005 

0.3356 

260.0000 
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TABLE  C6.   High-Intensity,  Vertical  Spatial  Correlation 
and  Structure  Function 


Thermistor  No, 


Variance  (°(T) 


1 

2 
3 
4 

5 
6 

7 

0.1800  X 

0.2462  x 
0.4578  X 
0.1231  x 

0.1226  X 

0.5298  X 
0.2655  x 

10' 
10" 
10" 
10' 
10" 
10" 
10" 

-5 

-3 
-3 
-1 
-2 
-2 
-1 

Correlation 
Function 

Structure 
Function 

0.0 

Lag 

1.0000 

0.0 

0.6971 

9.6513 

5.0000 

0.3813 

9.6565 

15.0000 

0.5125 

0.3485 

20.0000 

1.7985 

1.7245 

25.0000 

0.01*17 

0.1730 

25.0000 

1.1782 

12.0107 

30.0000 

0.8521 

2.9908 

45.0000 

8.3755 

9.2314 

55.0000 

3.7301 

15.2006 

80.0000 

2.9532 

25.7984 

85.0000 

1.9051 

18.2240 

100.0000 

0.2184 

0.1375 

135.0000 

0.2522 

9.7420 

150.0000 

0.2551 

0.6906 

155.0000 

0.0156 

0.3438 

160.0000 

0.0156 

10.0158 

175.0000 

0.0170 

0.9675 

180.0000 

0.3564 

3.8096 

180.0000 

0.0179 

4.2871 

205.0000 

1.0541 

19.7533 

235.0000 

0.0543 

21.5538 

260.0000 
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